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Abstract 

We coupled first principles calculations and the quasiharmonic approximation combined 
with experiments (X-Ray diffraction and dilatometry measurements) to determine the thermal 
properties of NiTiSn (half-Heusler) and Ni2TiSn (Heusler) compounds. These properties are 
important especially if they are to be used in thermoelectric applications. First, the calculation 
of their mode Gruneisen parameter shows that it is positive throughout the first Brillouin zone. 
This suggests that these compounds undergo a regular thermal expansion. Then, the calculation 
of the Ni2TiSn thermal expansion shows an excellent agreement, even in the high temperature 
range, with our high energy powder X-Ray diffraction measurements (ESRF) and dilatometry 
experiments. In the case of NiTiSn, this agreement is less impressive. This is due to stronger 
phonon-phonon interactions that are not considered within the quasiharmonic approximation, 
but also to the difficulty of making high-quality NiTiSn samples. Finally, the constant-pressure 
and constant-volume heat capacities have been calculated for both compounds and compared 
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with the experimental data reported in the literature. In particular, we have decomposed the 
constant-volume heat capacity of Ni2TiSn into a purely electronic and a phonon-mediated 
contribution, and we discuss each of them. 

Keywords: Thermal expansion, Griineisen parameters, heat eapaeity, Heusler eompounds, first 
principles, quasiharmonic approximation, X-Rays, dilatometry 

I. INTRODUCTION 

Heusler and half-Heusler compounds are the archetypes for multifunctional materials and their 
versatility increases constantly.® They generally crystallize into cubic structures (Fm3m for the 
Heusler and F43m for the half-Heusler are the usual space groups®), which increases their theo¬ 
retical attractivity. In addition, they are made of environmental friendly elements which permits to 
use them in practical applications. In this paper, our work is focused on the Ni 2 TiSn Heusler and 
NiTiSn half-Heusler compounds which have good thermoelectric properties, properties that are of 
interest to us. In previous studies, we have tackled through ab initio simulations the electronic 
properties® and the thermodynamic/mechanical® properties of these compounds since a compre¬ 
hensive study of these materials is still missing even though it is necessary if one wants to use them 
in real thermoelectric devices. 

The thermoelectric effect is a direct conversion of a temperature gradient into an electric current 
via the Seebeck effect. Thus, thermoelectric generators contribute to sustainability through scav¬ 
enging of waste heat or heat sources.® This implies that an accurate knowledge of the behavior 
of the thermoelectric materials as a function of temperature is necessary. Such a knowledge is 
missing for the Heusler type materials considered in our study. For the half-Heusler, Wee et a/.® 
have calculated the Griineisen parameters of the acoustic phonon modes to do a rough estimate 
of its lattice thermal conductivity, and Jung et a/.® have measured its thermal expansion up to its 
decomposition point. 

We propose here to couple the quasiharmonic approximation and the density functional per- 
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turbation theory to study the thermodynamie properties of NiTiSn and Ni 2 TiSn for temperatures 
up to 1500 K (elose to the melting temperature of NiTiSn). For this purpose, we first compute 
the mode Gruneisen parameter of the acoustic and optical phonon modes and the constant-volume 
heat capacity for both compounds. These quantities are indeed of special interest to clearly under¬ 
stand their thermal expansion mechanism. Then, we derive their constant-pressure heat capacity. 
In the case of Ni 2 TiSn, the latter is decomposed into a purely electronic and a phonon-mediated 
contribution, and we discuss each of these contributions. Finally, our calculations will be com¬ 
pared with experimental data. In this context, we measured the thermal expansion of Ni 2 TiSn for 
which no such data exist in the literature using temperature dependent powder X-Ray diffraction 
measurements from 80 to 673 K on laboratory equipment and at the European Synchrotron Re¬ 
search Facility (ESRF) as well as dilatometry experiments up to 873 K. Although the experimental 
thermal expansion of NiTiSn has been reported by Jung et al.P^ we revisited this work using new 
X-ray measurements for the sake of completeness and to compare our data with the ones already 
published. 

The paper is organized as follows: in section two, we present the experimental conditions 
whereas we detail the theoretical framework in section three. In section four, we show and discuss 
the results, while section five draws the major conclusions of our work. 

II. EXPERIMENTAL CONDITIONS 

The elements (Ti (99,7%, 90/rm), Ni (99,9% , 5 /im) and Sn (99,85% , 90 jim) were weighed in 
quantities corresponding to the desired stoichiometry. Eor the synthesis of Ni 2 TiSn, the elements 
were mixed and compacted at 375 MPa during 30 minutes under primary vacuum. Ingots of 
Ni 2 TiSn were then prepared by arc melting. The NiTiSn ingots prepared by arc melting as well 
were sealed into silica tubes and subsequently annealed at 1073K for ten days in order to improve 
their quality. The resultant materials were investigated by X-ray diffraction (XRD) using the Cu 
Kq;- radiation in the 26 range from 20 to 60°. The observation of the morphology of the materials 
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coupled to the quantitative chemical analysis was performed on a scanning electron microscope 
FEI Quanta 200, resolution 3 nm, vacuum, 30kV for the secondary electrons and 4 nm for the 
backscattered electrons, coupled with a microprobe EDS (Energy Dispersive X-ray Spectroscopy) 
Oxford Instrument XMax with a detector of 50 mm^. In each case, the desired phase was obtained. 

Concerning the Ni 2 TiSn material, measurements were made both at low and high temperatures. 
Eor low temperature experiments (performed at the European Synchrotron Radiation Eacility - 
ESRE, Erance), a fine powder of Ni 2 TiSn was loaded into a glass capillary of 0.5 mm diameter. 
The capillary was sealed and cooled from room temperature to 80 K and then heated at a rate 
of 30 K/hour up to 250 K, while synchrotron powder-diffraction data were collected in situ. X- 
Ray diffraction was performed on the Swiss-Norwegian Beam Eine (BMOIA) at the ESRE. All the 
diffraction patterns were collected using a monochromatic beam with the wavelength of 0.70814 A 
and a PIEATUS 2M detector. The sample-detector distance (343.71 mm) and parameters of the 
detector were calibrated using NIST standard EaB6. Two-dimensional diffraction images were 
integrated using the Eit2D software.^ The temperature was controlled with an Oxford Cryostream 
700-I-. During each collection time (120 s per image) the capillary was rotated by 60° in the same 
angular interval. X-ray measurements as a function of temperature in the temperature range 293 to 
673 K were also realized for Ni 2 TiSn (performed at the “Eaboratoire de Physique du Solide” (EPS) 
in Paris-Orsay) with a heating rate of 5 K/min. The X-ray set-up consists of a diffractometer in 
the normal-beam geometry. The sealed silica capillary containing the powder is fixed on a holder 
rotating perpendicularly to the incident beam. It is mounted on a rotating anode generator (copper 
anode: A = 1.542 A) equipped with a doubly curved graphite monochromator. The sample-holder 
was enclosed into a furnace. The diffraction pattern was recorded on an image plate detector 
(MAR 345) located at a distance of 425 mm from the sample and tilted of 20° from the incident 
beam direction. Einally we measured also the thermal expansion at the PRIME Verre company 
(Montpellier) with a ADAMEE-EHOMARGY model DI.10.2 dilatometer for rectangular-shaped 
samples cut from the ingots between room temperature and 873 K under air, with a heating rate of 
3 K/min. Erom the obtained thermal expansion curves, the average volumetric thermal expansion 
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coefficient was calculated with uncertainties of 0.1 MK~^. 

For the NiTiSn material, X-Ray analyses as a function of temperature were performed at the 
ICGM (Montpellier) using a temperature PANalytical X’Pert Pro Philips apparatus. The unit was 
equipped with an Anton Paar HTK 1200 high temperature chamber to work under vacuum (pres¬ 
sure V.IO^'^ mbar) and for temperatures between 293 and 673 K. The NiTiSn powders were first 
mixed with a alumina in an ethanol solution in order to obtain an homogeneous repartition of the 
powder in the sample holder, a alumina was used in these experiments as reference because the 
dilatation coefficient of this material as a function of temperature is well-known and served thus to 
correct the dilatation phenomena inherent to the apparatus. Acquisition was realized between 20 
and 90° (20) with a step size 0.0130° and a scanning time of 600 s. Dwell time was kept at 10 min 
for stabilizing the temperature before recording the X-rays. 

III. THEORETICAL FRAMEWORK 

A. Thermal expansion coefficient 

The equilibrium volume of a crystal at a given temperature T and in the absence of any applied 
pressure is obtained by minimizing the Helmholtz free energy with respect to all possible internal 
degrees of freedom: 

F{V, T) = EQ{V)+Fei{V, T) +F^,(y, T), (1) 

where V is the unit cell volume, Eq is the ground state (T =0 K) total energy of the crystal, Egi is 
the electronic free energy contribution, and Epi^ is the vibrational free energy which comes from 
the phonon contribution. In the quasiharmonic approximation (QHA), the latter term is written as 

follows 


Fph{V,T) 


^ 7>q 


1 — exp 


( ^to(y,q,y) Y 

V ksT J\^ 


( 2 ) 
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where the sums run over all allowed waveveetors q in the first Brillouin zone and over all phonon 
branehes j, kg is the Boltzmann eonstant, fi is the redueed Planek eonstant, and ( 0 ( 7 , q,y) is the 
frequeney of the phonon with waveveetor q in braneh j, evaluated at eonstant-volume V. The first 
and seeond terms of [^represent the zero-point and thermal energies of the phonons, respeetively. 
Fph ean be obtained by ealeulating the phonon dispersion relations and the eorresponding phonon 
density of states (DOS) from density funetional perturbation theory. The QHA assumes that the 
phonon frequeney is only a funetion of volume and is temperature-independent. Thus, the eleetron- 
phonon eoupling is negleeted. 

For the eleetronie eontribution to the free energy, Fgi = Egi — TSgi, the eleetronie energy due to 
thermal eleetronie exeitations is given by:^ 



( 3 ) 


where n(e, V), /(e), and £f represent the eleetronie density of states, the Fermi-Dirae distribution, 
and the Fermi energy, respeetively. The eleetronie entropy is formulated as: 


poo 

Sgi{VJ) = -kB / n(£,y)[/(e)/n/(e) + (l-/(£))/«(!-/(e))]Je 
Jo 


( 4 ) 


These eontributions are determined from the ealeulated eleetronie DOS. 

By definition, the volumetrie^ thermal expansion eoeffieient of a erystal is given by:^ 



( 5 ) 


where the subseript P implies that the temperature derivative is taken at eonstant pressure. For 
eubie symmetry, the minimization of[^with respeet to the volume leads to an alternative expression 
for av eomposed of the sum of a vibrational term (ct^*) and an eleetronie term 



( 7 ) 




where B is the bulk modulus and 7j(q) is the mode Griineisen parameter defined as: 


_ / aZn[m(7,q,y)] 

V dlnV yo 


( 8 ) 


Here, the subscript "0" indicates a quantity taken at the ground state lattice parameter. The vibra¬ 
tional specific heat at constant-volume is obtained using the calculated phonon DOS as follows: 




,/q 


V 2kfir 


( 9 ) 


and the electronic specific heat at constant-volume can be obtained from:^ 


rei -T I 
” “ ' dT 


( 10 ) 


The specific heat at constant-pressure can be determined by using the relation: 


Cp-C^^al{T).B{T).V{T).T, 


( 11 ) 


where Cy = Cy ^ -I- is the total specific heat at constant-volume. 


B. Computational details 

Density functional theory (DFT) based calculations were performed using the ABINIT packageff^ 
and the generalized gradient approximation (GGA) parametrized by Perdew, Burke and Emzer- 
hof (PBE).^The all-electron potentials were replaced by norm-conserving pseudopotentials gen¬ 
erated according to the Troullier and Martins scheme.^ Ni(3J^, 4s^), Ti(3J^, 4s^) and SnCSj'^, 
5p^)-electrons were considered as valence states in the construction of the pseudopotentials. The 
electronic wavefunctions were expanded in plane-waves up to a kinetic energy cutoff of 65 Ha. In- 
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tegrals over the Brillouin zone were approximated by sums over a 8 x 8 x 8 mesh of speeial ^-points 
aeeording to the Monkhorst and Paek seheme.^ A Fermi-Dirae seheme with a smearing width 
equal to 0.01 Ha was used for the metallic occupation of Ni 2 TiSn. Phonon dispersion curves were 
interpolated over a 4 x 4 x 4 (^r-points grid according to the scheme described by Gonze et 
whereas a denser 120 x 120 x 120 grid was employed for the calculation of the thermodynamic 
properties. 

IV. RESULTS AND DISCUSSION 

A. Structural properties 

In order to get the temperature dependence of the lattice parameters of NiTiSn and Ni 2 TiSn, we 
have calculated the total free energy at temperature points from 0 to 1500 K with a step of 15 K and 
for five volumes. At each temperature point, the equilibrium volume and the isothermal bulk mod¬ 
ulus, B{T) =V > is obtained minimizing the free energy from a Birch-Murnaghan equation 

of stateResults are shown inwith the experimental ones. Both for the calculation and the 
experiment, the equilibrium lattice parameters (top panel) show a regular thermal expansion: they 
increase with increasing temperature. The lattice parameters are sligthly overestimated by our cal¬ 
culations within 1-2% with respect to the experiments, as usual with GGA exchange-correlation 
functionals.^ The dependence of the isothermal bulk modulus as a function of temperature is dis¬ 
played in the bottom panel. For both compounds, it is almost a linearly decreasing function as the 
temperature increases. However, the slope of these lines is larger for the metallic compound than 
for the semiconducting one. As a consequence, these lines intersect near the melting temperature 
of NiTiSn (r,„ =1453 Below this temperature, the NiTiSn bulk modulus remains the lowest. 
This is due to the Ni-vacancies which make its structure softer to an hydrostatic pressure. To our 
knowledge, there is no experimental data to check our predictions. 
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B. Mode Griineisen parameters 

[^reports the mode Griineisen parameter of NiTiSn and Ni 2 TiSn as defined by along some high 
symmetry direetions. This ealeulation requires the knowledge of the phonon dispersion eurves at 
two additional unit cell volumes. They are derived from the equilibrium volume by straining it by 
±3%, and reoptimizing the atomic positions. The dispersion curves are discontinuous at the zone- 
center as a consequence of the polarization dependence of the sound velocities. Indeed, because of 
the vanishing of the acoustic frequencies at the F-point, the dispersions of the acoustic branches are 
discontinuous at the Brillouin zone center, and the value of the Griineisen parameters in the limit 
q —)■ 0 depends, for such modes, on the direction of q. We observe that the Griineisen parameters 
are positive throughout the Brillouin zone for all branches. This suggests from Eq. (7) that there 
is no anomalous negative thermal expansion at low temperatures in these compounds since the 
heat capacity and the bulk modulus are always positive. While no clear dominant mode Griineisen 
parameter is observed for NiTiSn, some branches have a high Griineisen parameter at the zone- 
center for Ni 2 TiSn. This contributes to the increase of the overall Griineisen parameter of Ni 2 TiSn 
at low temperature, since at these temperatures only the very low-frequency phonons contribute to 
the summation in[^ whereas high-frequency phonons decay exponentially. The overall Griineisen 
parameter of Ni 2 TiSn is 7=2.08, higher than the one of NiTiSn (7=1.60). The most positive 
Griineisen parameters in Ni 2 TiSn are located near the F-point and their frequencies are below 
150 cm^^ These modes mainly involve motions of Ni-atoms (see for instance the case of the T 2 g 
Raman mode calculated at 122 cm ^ in Ref.^. 

C. Thermal expansion coefficients 

In 1^ we show the temperature dependence of the volumetric thermal expansion of NiTiSn and 
Ni 2 TiSn. This thermal expansion has been calculated using Eq. (7) with temperature-dependent 
values of the bulk modulus and the unit cell volume obtained after minimization of the free energy 
from a Birch-Murnaghan equation of state. The experimental thermal expansion was obtained from 
Eq. (5) with a derivative evaluated numerically using three consecutive points along the curve giv- 
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ing the experimental lattiee parameter as a funetion of temperature. This methodology is expeeted 
to be more aeeurate (but less smooth) than a proeedure involving a polynomial fit of this eurve and 
then making the derivation, beeause it is independent of the degree of the polynomial used for the 
fit. 

The deeomposition of the thermal expansion into a vibrational and an eleetronie eontribution is 
also reported on this figure for Ni 2 TiSn. We observe that the vibrational eontribution mainly dom¬ 
inates the thermal expansion while the inelusion of the eleetronie eontribution leads to a minor 
eorreetion up to 1500 K. In the ease of NiTiSn (semieondueting eompound), eleetronie eontribu- 
tions are negligible at the first order sinee the eonsidered temperatures are signifieantly lower than 
the ealeulated energy band gapS!(£’^ = 0.49 eV). 

NiTiSn has a lower thermal expansion eoeffieient than Ni 2 TiSn in the whole range of tempera¬ 
ture due to its smaller Griineisen parameter. This is also a eonsequenee of the low phonon density 
observed at low frequeneies (below 100 em^^) in its vibrational DOS with respeet to Ni 2 TiSn.® 
It is important to bear in mind that the QHA used in the ealeulations of ay negleets anharmonie 
effeets, arising from phonon-phonon interaetions, whieh means that it ean break down at high 
temperatures where sueh effeets ean beeome signifieant. For Ni 2 TiSn, the agreement is very good 
between the ealeulations and the experiments in the reeorded temperature range (80-873 K) even if 
the high temperature XRD measurements are slightly overestimated. The good agreement between 
experiments and simulations suggests that the QHA remains valid above room temperature and at 
least up to 1000 K in this eompound. The agreement is less impressive between our ealeulated 
and high temperature experimental values for NiTiSn even if the ealeulated eurve goes through 
the experimental points. A possible explanation is that the phonon-phonon interaetions (normal 
and umklapp proeesses) are more important in this eompound than in Ni 2 TiSn. However, these 
interaetions alone ean not explain this lesser agreement. Indeed the NiTiSn sample is probably 
not monophasie beeause it is very diffieult experimentally to make high-quality samples, while 
the synthesis of Ni 2 TiSn is relatively easy and straightforward. Thus, phonon-defeet interaetions 
should probably also be taken into aeeount for NiTiSn. 
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Concerning particular numbers, at room temperature our calculations and dilatometry experi¬ 
ments give: ay — +40 for Ni 2 TiSn. For NiTiSn, we find ay — +30 in the calculations 
while the experimental value (XRD) is slightly higher: ~ +31 At 750 K, the temper¬ 

ature at which the thermoelectric efficiency of NiTiSn is maximum,!^ we obtain a value of ay ~ 
+33 MK ^ close to the thermal expansion coefficient of iron at room temperature. The main errors 
in these calculated values of the thermal expansion are mainly related to: (i) the calculated equi¬ 
librium lattice constant that is slightly different from the experimental measurements (see[^, (ii) 
the type of exchange-correlation functionals used in DFT, and (iii) the conditions where the QHA 
remains valid above room temperature. 

Note that in the literature, the average linear thermal expansion can also be defined by:l^ 


(«) 


ave 


1 L(T)-Lrt 
Lrt T — Trt 


( 12 ) 


where Lrt is the sample length at the room temperature reference and this definition was used to 
determine experimentally ay of NiTiSn by Jung et a/.^^With this new definition, we have calculated 
the thermal expansion of NiTiSn using Eq. (7) with the values of the bulk modulus and unit cell 
volume fixed at 298 K. In our calculations, these values are B{29SK) = 114 GPa and V{29SK) = 
54.6 A^. In the inset of[^we compare the so-obtained calculated values with the experimental 
results of Jung et al.^ and with our own experimental points determined with this equation. We 
observe an overall good agreement between the three sets of experimental data and the calculated 
curve. Using the same methodology for Ni 2 TiSn for the calculated and experimental points, we 
found that the thermal expansion of both compounds grows rapidly up to 350 K and then becomes 
nearly constant above this temperature. Thus, the main differences between the thermal expansion 
calculated using Eq. (7) and Eq. (12) appear roughly above room temperature (see|^. 


11 



D. Heat capacities 


Once the phonon speetrum over the whole Brillouin zone is available, the vibrational heat eapaeity 
at eonstant-volume (Cv^) ean be ealeulated by while the eleetronie eontribution to the heat 


eapaeity at eonstant-volume (C®') ean be obtained from the eleetronie DOS by using 10 The 
temperature dependenee of the speeifie heats at eonstant-pressure Cp of both eompounds have been 
ealeulated using The different heat eapaeities and eontributions are plotted in|^in the ease of 
Ni 2 TiSn, whereas only the vibrational eontribution to Cv is eonsidered for NiTiSn in First, 
the Cv^ eontributions tend to the elassieal Dulong and Petit eonstant as the temperature inereases. 


Cv^ = liNkB where N is the number of atoms, while Cf and Cp still inerease. In partieular, Cv'\T —)■ 
oo) =99.77 J.mol"for Ni 2 TiSn and cf{T oo) =74.83 J.moDfor NiTiSn. Then, for 
the thermal eleetronie eontributions to the speeifie heat of Ni 2 TiSn, we find that is not negligible 
at high temperatures though signifieantly smaller than Cy^. This eharaeter ean be understood 
from the high eleetronie DOS observed near the Fermi level.^ For NiTiSn, the small differenee 
between Cp and Cy at high temperatures is a eonsequenee of its smaller thermal expansion. In 
eontrast to the high temperature range, the differenee between Cp and Cy ean be negleeted for both 
eompounds at low temperature (below ~ 300 K). For Ni 2 TiSn, we observe for Cp a remarkable 
agreement between our ealeulated values and the experimental^ ones above room temperature. 
This highlights the reliability of the ealeulation of the thermal expansion and the validity of the 
QHA above room temperature in this eompound. Below this temperature, only measurements of 
the speeifie heat below 25 K have been reported in the literature to our knowledgej ^^ ^ ^^ l leaving 
the 25-300 K range unexplored. In the ease of NiTiSn, this agreement is also satisfaetory with the 
experimental data obtained by Zhong^^ up to 400 K. There is no experimental values reported in 
the literature above room temperature. 


-.ph. 
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V. CONCLUSIONS 


We have determined via first prineiples ealeulations and experimental measurements ineluding 
three different methods (high energy X-Rays (ESRF), standard X-Rays (Montpellier) and dilatom- 
etry) the thermal eharaeteristies of NiTiSn and Ni 2 TiSn. First, the ealeulated dispersion eurves of 
the mode Gruneisen parameter are positive throughout the Brillouin zone, suggesting that these 
eompounds undergo a regular thermal expansion. No elear dominant mode Gruneisen parame¬ 
ter is observed for NiTiSn whereas some branehes have a high Gruneisen parameter at the zone 
eenter for Ni 2 TiSn. As a eonsequence, the overall Gruneisen parameter of Ni 2 TiSn is the highest. 
Then, the ealeulated thermal expansion eoeffleients are in exeellent agreement with the experimen¬ 
tal values in the reeorded temperature range (80-873 K) for Ni 2 TiSn. In the ease of NiTiSn, this 
agreement is less impressive at high temperatures. These deviations ean be attributed to stronger 
phonon-phonon interactions that are neglected within the QHA, but also to phonon-defect interac¬ 
tions due to the low quality of the sample. The constant-volume heat capacity has been calculated 
for both compounds and compared with the experimental data reported in the literature. In partic¬ 
ular, we decomposed that of Ni 2 TiSn into a purely electronic and a phonon-mediated contribution, 
and we discussed each of them. Finally, we derived their constant-pressure heat capacities. The re¬ 
markable agreement found for Ni 2 TiSn with the experimental data highlights the reliability of the 
calculation of its thermal expansion and the validity of the QHA above room temperature in this 
compound. This exhaustive study of the thermal properties of NiTiSn and Ni 2 TiSn is the first step 
in the determination via ab initio methods of the thermal conductivity of these Heusler materials. 
Such completely first principles methods do not exist yet even though obtaining reliable thermal 
conductivities is fundamental to predict new and efficient thermoelectric materials. 
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Figure 1: Temperature dependence of the lattice parameter (top) and the bulk modulus (bottom) 
for Ni 2 TiSn and NiTiSn. Full squares represent the experimental data from Jung et al^ 
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Figure 2: Calculated dispersion curves of the mode Griineisen parameter (7;(q)) for Ni 2 TiSn (left) 
and NiTiSn (right) along some high symmetry lines. 
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Figure 3: Temperature dependence of the volumetric thermal expansion of NiTiSn and Ni 2 TiSn 
obtained using Eq. (5): comparison between XRD, dilatometry and calculations. In the inset of 
the left panel we compare our XRD and calculated values of NiTiSn obtained using Eq. (12) to 
the experimental data of Jung et also reported in their paper using Eq. (12). Note: IMK^^ is 
equivalent to 1 x 10 ^ 
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Figure 4: Temperature dependence of the volumetric thermal expansion of NiTiSn and Ni 2 TiSn 
calculated following Eq. (7) (theoretical determination) and Eq. (12) (experimental determination). 
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Figure 5: Temperature dependence of the Ni 2 TiSn specific heat at constant-volume (Cy) and 
constant-pressure {Cp). Electronic {Cf) and phonon (Cf^) contributions to the specific heat at 
constant-volume are also reported. The open circles are the experimental values from Nash.^^ 
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Figure 6: Constant volume (Cy^, vibrational part), and constant pressure (Cp) specific heat of 
NiTiSn versus temperature. The open circles are the experimental values from Zhong 
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